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Abstract: 

We investigate the possibility that gamma-ray bursts are powered by a central 
engine consisting of a black hole with an external magnetic field anchored 
>• I in a surrounding disk or torus. The energy source is then the rotation of 

the black hole, and it is extracted electromagnetically via a Poynting fiux, a 
CSJ ' mechanism first proposed by Blandford and Znajek (1997) for AGN. 

p^ ! Our reanalysis of the strength of the Blandford-Znajek power shows that 

Q\ \ the energy extraction rate of the black hole has been underestimated by a 

^ I factor ten in previous works. Accounting both for the maximum rotation 

1^ • energy of the hole and for the efficiency of electromagnetic extraction, we 

find that a maximum of 9% of the rest mass of the hole can be converted to 
a Poynting fiow, i.e. the energy available to produce a gamma-ray burst is 
^1 1.6 X 1O^^(M/M0) erg for a black hole of mass M. We show that the black 

^ ' holes formed in a variety of gamma-ray burst scenarios probably contain the 

required high angular momentum. 
^ ■ To extract the energy from a black hole in the required time of ^ 1000 

^ . s a field of 10^^ G near the black hole is needed. We give an example of 

a disk-plus-field structure that both delivers the required field and makes 
the Poynting fiux from the hole dominate that of the disk. Thereby we 
demonstrate that the Poynting energy extracted need not be dominated by 
the disk, nor is limited to the binding energy of the disk. This means that 
the Blandford-Znajek mechanism remains a very good candidate for powering 
gamma-ray bursts. 
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1 Introduction 

Gamma ray bursts presently provide great excitement in astronomy and as- 
trophysics as optical observations by way of many instruments give consider- 
able detail of the history of each burst. We are concerned here with the prodi- 
gious energy in each burst, the estimate for GRB 971214 being ^ 3 x lO^^ergs 
[|5], although this could be diminished if considerable beaming is involved in 
the central engine, as we will discuss. 

Amazingly, 2 x 10^'^ergs is just the rest mass energy of our sun, so it 
seems immediately clear that the central engine for the GRB must be able 
to extract a substantial fraction of the rest mass energy of a compact object, 
neutron star or black hole, and convert it into energy of GRB. 

The second criterion for the central engine is that it must be able to 
deliver power over a long interval up to ~ 1000 seconds, since some GRB's 
last that long although other GRB's last only a fraction of a second. It must 
also be able to account for the vast diversity in pulses, etc., or, alternatively, 
one must have a number of diverse mechanisms. 

We believe the need to deliver power over the long time found in some 
bursters to be the most difficult requirement to fulfill, since the final merger 
time of the compact objects is only a fraction of a second and it is difficult 
to produce a high energy source of, e.g., z/i/-collisions that goes on for more 
than two or three seconds. 

For many years mergers of binary neutron stars were considered to be 
likely sources for the GRB's. The estimated merger rate in our Galaxy 
of a few GEM ^ is of the right order for the occurence of GRB's. The 
possible problem with binary mergers might be the ejected materials during 
the merging processes. Not more than ~ 1O~^M0 of baryons can be involved 
in the GRB, since it would not be possible to accelerate a higher mass of 
nucleons up to the Lorentz factors F ~ 100 needed with the energies available. 

We find the merger of a neutron star with a black hole to be a particu- 
larly attractive mechanism. The baryon number "pollution" problem can be 
solved by the main part of the baryons going over the event horizon. In the 
Blandford-Znajek mechanism[0] we wish to invoke, a substantial proportion 
of the rotational energy of the black hole, which will be sent into rapid ro- 
tation by swallowing up the neutron star matter, can be extracted through 



1 GEM = 1 Galactic Event per Megayear 



the Poynting vector. The rate of extraction is proportional to the square of 
magnetic field strength, B^, as we shall discuss, so that power can be fur- 
nished over varying times, depending upon the value of B. With substantial 
beaming, we estimate that B ~ lO^^G would be sufficient to power the most 
energetic GRB's with ~ lO^^ergs. 

Recently at least three magnetars, neutron stars with fields B ~ 10^^ — 
lO^^G, have been observed. Their visible lifetime is only a few thousand years, 
because neutron stars slow down by emission of magnetic dipole radiation 
and join the "graveyard" where they no longer emit pulses. The time of 
observability is proportional to B~^, so the number of these high magnetic 
field stars may not be an order of magnitude less than the garden variety 
lO^^G neutron stars. It is also possible that existing magnetic fields can be 
increased by the dynamo effect. 

Failed supernovae were suggested by Woosley(1993)|Q as a source of 
GRB. In this case the black hole would be formed in the center of a massive 
star, and surrounding baryonic matter would accrete into it, spinning it up. 
This mechanism is often discussed under the title of hypernovae[^]. 

More recently Bethe and Brown(1998)|P] found that in binary neutron 
star evolutions, an order of magnitude more low-mass black-hole, neutron- 
star binaries were formed than binary neutron stars. The low-mass black-hole 
mass of ~ 2AMq looks favorable for the Blandford-Znajek mechanism. 

In some calculations which begin with a neutron star binary, one of the 
neutron stars evolves into a black hole in the process of accretion, and the 
resulting binary might also be a good candidate for GRB's. In any case, 
there are various possibilities furnished by black-hole, neutron-star binaries. 

In this paper we will discuss the Blandford-Znajek mechanism in quan- 
titative detail. In section 2, an overview of the proposed central engine for 
gamma ray bursts using Blandford-Znajek process is given. The Blandford- 
Znajek process and the evolution of the black hole will be discussed in detail 
in section 3 and section 4 respectively. The structure of ambient magnetic 
field surrounding black hole and accretion disk is discussed in section 5. In 
section 6, we will give a rough estimation of the possible angular momentum 
of the black hole which might emerge as a final compact object during the 
merging or collapsing processes. The possible constraint from the surround- 
ing accretion disk is discussed in section 7 and the results are summarized 
and discussed in section 8. 



2 Overview of the Proposed Central Engine 
for GRB 

Two decades ago Blandford and Znajekp| ^ proposed a process(BZ) in which 
rotational energy of black hole can be efficiently extracted. If there are suf- 
ficient charge distributions around the black hole to provide the force-free 
condition, then the magnetic field lines exert no force and corotate rigidly 
with the rotating black hole. The induced current loops loops which pass 
along the black hole's stretched horizon feel the forces by the magnetic field 
supported by the environment. Hence these forces give magnetic braking 
of the black hole rotation. The maximum amount of energy which can be 
extracted out of the black hole without violating the second law of thermo- 
dynamics is the rotational energy, Ej-ot, which is defined as 

Erot = Mc^-MirrC^ (l) 

(2) 
where 
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where Ah and Sh are surface area and entropy of a black hole respectively^ 
and kB is Boltzman's constant. The rotational energy of a black hole with 
angular momentum J is a fraction of the black hole mass M, 

Erot = f{a)M^ (5) 

/(a) = l-^i[l + v/r^ (6) 

where a = -^§q is the rotation parameter. For a maximally rotating black 
hole (a = 1), / = 0.29. In BZ, the efficiency of extracting the rotational 



■^For a solar mass black hole, Ah ~ lO^^cm^, S'/f ~ lO^^/cs. The Planck mass is 



energy is determined by the ratio between the angular velocities of black 
hole, Qh, and magnetic field angular velocity flp, 

The rest of the rotational energy is dissipated into the black hole increasing 
the entropy or equivalently irreducible mass^ The total BZ energy available 
is 

Ebz = 1.8 xlO'^enf (a) i^)erg (8) 

For the optimal processes e^ ~ 0.5 ||^. 

Since the energy transport is in the form of Poynting fiow in BZ the out- 
going energy fiux from the black hole is basically B^c. An order of magnitude 
calculation is now in order. There are basically three parameters: Mass of 
the black hole M and magnetic field on the horizon B, which are dimension- 
ful, and angular momentum parameter of the black hole a. The time scale 
for the BZ process can be calculated by the ratio of the black hole mass to 
the output power from the black hole surface B^R^c 
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B^R^c B^M^G^ 

2.7 X 10^(Hi!^^!^)^(^). (9) 
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where Mq^ {10^^ gauss)'"^-^ = 2.7 x lO'^s and the radius of horizon, R, is 
taken to be ~ ^^. Also the outgoing Poynting power is 

Psz - B^R'c = 6.7 X 10^°( ,n/ ?i^?^rg/s. (10) 

The fiuence of the recently observed GRB971214[|l) corresponds to E^ = 
lO^^-^{^)erg which is consistent with Ebz- That of GRB990123 may be 
as large as E^ = 3.4 x 10^^(^)|^. This suggests that if a strong enough 
magnetic field (~ 10^^) on the black hole can be supported by the surrounding 
material(accretion disk) the BZ process is a good candidate to provide the 

^The mass of a nonrotating black hole itself is its irreducible mass 



powerful energy of the GRB in the observed time interval up to 1000s , which 
is comparable to the BZ time scale tbz- 

In recent years a black hole plus debris torus system(or accretion disk) 
has been considered to be a plausible structure for the GRB central engine ||]. 
The presence of the accretion disk is important for the BZ process because 
it is the supporting system of the strong magnetic field on the black hole, 
which would disperse without the pressure from the fields anchored in the 
accretion disk. Recent numerical calculations show that accretion disks 
formed by various merging processes are found to have large enough pressure 
such that they can support ~ lO^^G assuming a value of the disk viscosity 
parameter a ~ 0.1, where a is the usual parameter in scaling the viscosity. 
This gives a relevant order of magnitude magnetic field on the black hole, 
which is, however, not considered to be much larger than magnetic field of 
the inner accretion disk |Tn][|n[]. The discovery that soft gamma ray bursts 



are magnetars [0 also supports the presence of the strong magnetic fields of 
~ lO^^G in nature. The identification by now of three soft gamma repeaters 
as strong-field pulsars indicates that there may be a large population of such 
objects: since the pulsar spindown times scale as B~^, we would expect to 
observe only 1 magnetar for every 1000 normal pulsars if they were formed at 
the same rate, and if selection effects were the same for the two populations. 
We see 3 magnetars and about 700 normal pulsars, but since they are found in 
very different ways the selection effects are hard to quantify. It is nonetheless 
clear that magnetars may be formed in our Galaxy at a rate not very different 
from that of normal pulsars. 

The life time of the accretion disk is also very important for the GRB time 
scale because it supports the magnetic field on the black hole. According to 
numerical simulations of merging systems which evolve eventually into black 
hole - accretion disk configuration the viscous hfe times are 0.1 - 150 s, 
which are not inconsistent with the GRB time scale. Also it has been pointed 
out[^ that a residual cold disk of ~ 1O~^M0 can support lO^^G, even after 
the major part of the accretion disk has been drained into the black hole or 
dispersed away. 

Recent hydrodynamic simulations of merging neutron stars and black 



holes|]13[ show that along the rotation axis of the black hole an almost baryon- 
free funnel is possible. This can be easily understood since the material above 
the hole axis has not much angular momentum so that it can be drained 
quickly, leaving a baryon-free funnel. Hence relativistically expanding jets 
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along the funnel, fueled by Poynting outflow which is collimated along the 
rotation axis[14], can give rise to gamma ray bursts effectively. It has been 
observed that the BZ process is also possible from the disk since the magnetic 
field on the disk is not much less than that on the black hole[|10[. However 



the energy outflow from the disk is mostly directed vertically from the disk 
where the baryon loading is supposed to be relatively high enough to keep 
the baryons from being highly relativistic. Therefore the BZ process from 
the disk can be considered to have not much to do with gamma ray burst 
phenomena. However the BZ from the disk could power an outflow with 
lower r, but nonetheless high energy, which could cause an afterglow at 
large angles. That would lead to more afterglows being visible than GRBs, 
because the afterglows are less beamed. The BZ mechanism can also play 
a very important role in the disk accretion because it carries out angular 
momentum from the disk. We will discuss this in more detail later. 

The structure we are proposing as a central engine of the GRB is a system 
of black hole - accretion disk(debris of the torus): 

The rotating black hole is threaded by a strong magnetic field. Along the 
baryon-free funnel relativistic jets fueled by Poynting outflow give rise to 
the GRB. The interaction between disk and black hole is characterized by 
accretion and magnetic coupling. We consider the BZ process only after the 
main accretion process is completed, leaving an accretion disk of cold residual 
material, which can support a strong enough magnetic field. 

3 Blandford-Znajek Process 

Consider a half hemisphere (radius R) rotating with angular velocity Q and 
a circle on the surface at fixed 6{ in the spherical polar coordinate system) 
across which a surface current I flows down from the pole. When the external 
magnetic field B is imposed to thread the surface outward normally, the 
surface current feels a force and the torque due to the Lorentz force exerted 
by the annular ring of width Rd6 is 

AT = -RsinexIRABO (11) 

= -J-BAAann. (12) 

= -^A* (13) 



where A\l/ is the magnetic flux through annular ring extended by dO with 
surface area Aann- We consider an axially symmetric situation. From this 
magnetic braking, we can calculate the rotational energy loss rate 

AProt = nxAT (14) 

= f!^A*. (15) 

Blandford and Znajekpj demonstrated that such a magnetic braking is 
possible, provided that the external charge distribution can support an elec- 
tric current with the magnetic field threading the horizon. The original for- 
mulation of Bland ford and ZanjekP] is summerized in Appendix 3. Macdon- 



ald and Thorne fTq] reformulated the Blandford-Znajek process using a 3+1 
dimensional formalism, in which the complicated physics beyond the hori- 
zon can be expressed in terms of physical quantities defined on the stretched 
horizon]^. It can be shown that the rotational energy loss due to magnetic 
braking can be obtained from eq.fpTSD by simply replacing Q and B with 

n^riH, B -^ Bh (16) 
where H denotes the quantity on the stretched horizon 0: 

AProt = fi/fxAT (17) 

= Qh^A^. (18) 

The current / induced by the black hole rotation and the angular velocity 
Qf{0) of the rigidly rotating magnetic field which is dragged by the rotating 
black hole can be determined together with the magnetic field by solving 
Maxwell's equations: 

F't;; = AttJ'' (19) 

with the force free- condition 

F'^V^ = 0, (20) 

where J'^ is a current density vector. The detailed structure of the magnetic 
field will be discussed in section 5. In the BZ process, however, the power 
which can be carried out as Poynting outflow along the field lines is 

APmag = ^F" A* (21) 

Ztt 



Then the rest of the rotational energy is used to increase the entropy(or 
equivalently irreducible mass) of the black hole. Therefore the efficiency of 
the BZ process which is defined as the ratio of Pmag to Prot is 

The ideal efficiency, e^ = 1, is meaningless because for Vtp = ^h the Poynting 
outflow itself is zero as can be seen below in eg. (P^) . The optimal power can 
be obtained at e^ = 1/2||15|. Then the rest of the rotational energy is used 
to increase the entropy(or equivalently irreducible mass) of the black hole. 
Now consider the loading region far from the black hole, onto which magnetic 
flelds out of the black hole anchor. In most of the cases we are interested 
in, the inertia of the loading region can be considered to be so large that the 
transported angular momentum cannot give rise to any substantial increase 
of the angular velocity of the loading region. Therefore the angular velocity 
of the loading region can be assumed to be zero and the power delivered by 
the torque, eq.(0), along the fleld line is given 

APl = npAT (23) 

= Qf^A^Ii (24) 

which is identifled as the BZ power for GRB: 

APbZ = AP^mag = APl (25) 

Since the current / induced by black hole rotation is given by[^] 

iie) = ^^inH~np{e))Co''BH (26) 

the rotational energy loss and the BZ power can be given by 

AP.„, = ^Mk^M^^fi^Avl/ (27) 

An 

APsz = ^H^H-^i^) -.^^^^^ (28) 

Aire 

where u) is a kind of cylindrical radius deflned in Boyer-Lindquist coordinates 
(Appendix 1). 
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To get the total power, eq.(p8D should be integrated from 6* ~ to Obz, 
up to which the magnetic field lines from the black hole anchored on to the 
loading region for GRB. As a first approximation, we put 

Obz ~ vr/2 (29) 

Buie) ^ <Bh> (30) 

and we get an optimal power 

P..^1.7xlO»a'^,iL).,_5^)V(/.)e.«/. (31) 

The details are given in Appendix 2. The rest of the rotational energy given 
by 

Ph = Prot — Pbz (32) 

= %^Pbz (33) 

= Pbz, for optimal case (34) 

is dissipated into the black hole increasing the irreducible part of the black 
hole mass(equivalently increasing the entropy of the black hole). 

For 6bz < < 7r/2, the magnetic field lines from the black hole can 
be anchored onto elsewhere than the loading region, for example, onto the 
inner accretion disk. Then the angular velocity of the field line is determined 
by the angular velocity of the disk VLd-, ^f = ^d- The innermost radius of 



the disk can be considered to be the marginally stable orbit, r^s |16]. For 



an orbit rotating in the same direction as the black hole rotation, we have: 

G 



c 



2 



M[3 + Z2- V(3-^i)(3 + Zi + 2Z2)], (35) 



Z, = l + {l-aY^[(l + ayr^ + (l-ay/% (36) 



Z2 = ^J3a^ + Zl (37) 

Here r^s becomes GM/c^ as a ^ 1 (extreme rotation). In this limit, the 
angular velocity of disk Qo ~ ^h- Then one can expect Qp ~ ^h so that 
there is no Poynting fiow to the disk. However in this case there is no BZ 
Poynting outfiow. For finite a < 1, we can have either Qh > ^d ot 
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VIh < ^D, which can be determined by solving eg. (|65D . However if we can 
assume that the power from/to the disk can be much suppressed compared 
to that of the loading region, we can ignore the magnetic coupling between 
the black hole and the accretion disk. For example, suppose the portion of 
black hole - threading magnetic fields anchored onto the disk is somehow 
suppressed, 9bz ~ 7r/2, then we can assume there is no significant energy 
and angular momentum feedback into the disk due to magnetic coupling. If 
it is not so we can discuss only the limiting case. 

4 Evolution of a black hole via the Blandford- 
Znajek process 

While the black hole is slowed down and a part of the rotational energy is 
carried out as Poynting outflow, the rest of the rotational energy increases 
the entropy of the black hole or its irreducible mass. The increasing rate of 
the irreducible mass is given by using eq.p7|) and 



^ = Prot-PL (38) 

= [ ^^""'^^^' co'EnA^. (39) 

J Atic 

The irreducible mass eventually becomes the mass of Schwarzshild black hole 
when it stops rotating. The difference between the initial Kerr black hole 
mass. Mo, and the final Schwarzshild mass, M, is the energy output from 
black hole via the Blandford-Znajek process. The evolution of a Kerr black 
hole is determined by the evolution of its mass and angular momentum given 
by 

Using eq.(p3[) and (^TD in the Blandford-Znajek process, the evolution of the 
mass and the angular momentum are related by 

dM ^ dJ 

= Qp (42) 

dt dt ^ ^ 
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For the optimal case, e^ = 1/2(^2^7' = VLh/2)^ we can obtain an analytic 
expression for the mass in terms of the angular momentum. With the angular 
velocity of a black hole expressed in the angular momentum given by in the 
natural units G = c = 1 

nn = f (43) 

2M3(1 + ^1 - jyM^) 



we get 



dM J dJ 



dt 


4M3(1 + ^1 - J^/M^ 


^dM' 


1 dJ^ 


dt 


l + ^l-.P/M^ dt 


Intoducing a new variable H defined as 


H -- 


^i-^/^' 


eq.(|45|) can be written 






^dM^ 1 dJ^ 
dt H dt' 


Using the identities 




d j^ 

dVM^' 


^ ^(-< 


dJ^ 
dt 


^ ir--(-^ 


J^/M^ 


= 2H-H^ 



4^ dt 



dt 



eq.(|47|) can be written 

1 dM _ (l-H) dH 
mUT ~ 2H^ dt ' 



which can be integrated analytically to give [p^] 



(44) 
(45) 



(46) 



(47) 



(48) 
H)'4 (49) 



(50) 



(51) 



M = Moe^^J^, (52) 

V H 
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where Ho and Mq represent the initial angular momentum and mass of the 



black hole. From eq.(46) one can see that H = 1 for the maximally rotating 
(a = 1) black hole and H = 2 for the nonrotating one. 

Consider the black hole initially maximally rotating, which is slowed down by 
the Blandford-Znajek mechanism in the optimal mode(f2^ = VLh/2Y . The 
final black hole mass is given by 

gl/4 

M = -^Mo (53) 

= OMMo. (54) 

We see that 9 % of the initial mass or 31 % of the rotational energy can be 
taken out to power the gamma ray burst from the maximally rotating black 
hole. The extracted energy is therefore less than a half of the initial rotational 
energy. This can be easily understood by noticing that the fraction of the 
rotational energy drops faster because of the decreasing total mass and at the 
same time increasing irreducible mass in eq.(l). For a = 0.5, M = 0.98Mo 
or 2 % of the initial mass can be used to power the gamma ray burst. The 
time dependence of the power can be obtained from eq.( p47| ) using eq.(^): 



dt 4 (^H-H^)M^Bj,-, (55) 

where f{h) is defined in Appendix 2 and H = yt/^- The initial rate for 
the maximally rotating black hole can be obtained by taking H = 1 and 
/ = vr — 2. Eq. (|55D can be written as 

1 dM TT - 2 , G^ 

-Bh-r. (56) 



M"^ dt 4 ^ c5 

Assuming there is no change in magnetic field which is supported by the 
environment, we get 

If we extrapolate eq.(|57D until the black hole stops rotating, M — > 0.91 Mq, 
we get the time scale of Blandford-Znajek process as 

'■''"^' (58) 



M,Bf,G' 



^The more general discussion has been given by I. Okamoto||T^ where C = l/e^ 
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Since there is no considerable change of the weighting factor j{K) in eq.(|55|) 
from the maximally rotating black hole to the nonrotating black hole(/ = 
2/3), Eq.(|59|) can be considered as a reasonable estimate of a time scale, 
which is consistent with the rough estimate given in eq.(P). 



5 Magnetic field and Force-free plasma 

The rotating black hole immersed in the magnetic field induces an electric 
field around the black hole|0. The electromagnetic field in the vicinity of 
rotating black hole immersed in the uniform magnetic field in free space has 
been obtained [0 [jT9| as a solution of the source-free Maxwell equation 



F'^:: = 0. (60) 

From the analytic expression (Appendix 2) which gives an asymptotically 
uniform magnetic field at infinity(r —* oo) 

B = Bz, r^oo, (61) 

one can see that the radial component of the magnetic field on the horizon, 

B^ = , 2 ^ f m/ .n,M - {a/cf){rl - {a/cfcosH) 
[rfj + [a/cY cos"^ oy 

+2{a/cfrH{rH - M) (1 + cos^ 9)] (62) 

(63) 

vanishes as the rotation of black hole approaches extreme rotation 

a -^ GM/c, th -^ GM/c^. (64) 

This is what has been observed [|l^ |^ as the absence of the magnetic flux 
across the maximally rotating black hole. One can also see that there is 
no outward Pointying flow, which means no outflow of energy to infinity. 
Essentially it is because of the abscence of the toroidal component of the 
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magnetic field due to the vacuum environment which is charge and current 
free space. In other words, there is no current on the stretched horizon on 
which the external magnetic field exerts torques to slow the black hole down 
so as to extract its rotational energy. Therefore it is neccessary to have a 
magnetosphere with charges and currents to extract the rotational energy of 
the black hole. The force-free magnetosphere around a rotating black hole 
has been proposed 0: 

F'f;: = AtxJ^" (65) 

where J^ is a current density vector and 

F'^'J^ = (66) 



which is the force- free condition. From eq. (|66D , we get the degenerate con- 
dition 

E-B = (67) 

In the case of a rotating black hole in charge-free space one can see using 
eqs. in Appendix 2 

E^ -B^ ^0 (68) 

It should be remarked that since the force-free condition is essential for the 
BZ process the expulsion of magnetic field on the rotating black hole demon- 
strated for the charge-free space cannot be directly addressed in the BZ 



process pod . 



To maintain the current flows charged particles (electrons and positrons) 
should be supplied by a source outside the horizon since the black hole it- 
self cannot provide the outgoing particle from inside the horizon. Blanford 
and Znajekpl proposed that the strong electric field induced by the rotating 
black hole can give rise to a spark gap in which sufficient charged parti- 



cles are created to supply the currents||21||. Another mechanism provides 
charged particles around black hole: electron-positron pair creation by neu- 
trino annihilation!^. Recent numerical studies of merging binary systems 



[p|][|23| which result in black hole-accretion disk configurations show that the 
power of electron-positron pairs by the annihilation of neutrinos and antineu- 
trinos which are radiated out of the disk is 

Kp ~ lO^^erg/s. (69) 
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This power is being poured into the space above the black hole for 0.01 — Is. 



If we divide it by the average neutrino energy < e^ >~ IQMeV [^, we get 
a rough estimate of the numbers of e^e~ pairs 

Npair ~ lO^Vs (70) 

or the charge producing rate for e^(or e~;the sum is always zero) is 

Qe = eNpan^ ~ W'^'C/s (71) 

The magnitude of the currents involved in the optimal BZ process can be 



obtained from eq. (|250|) 



/ r. 12!^ ^ 10-C/s (72) 

V ^H 

for a black hole with solar mass threaded by a lO^^G magnetic field. Rh 
is the surface resistance of the horizon, 377 Ohm. From the comparison 
of eq.(^) with (^), the neutrino annhilation process produces orders of 
magnitude more than enough pairs to keep the necessary currents for the 
optimal BZ process. The possible effects of neutrino annihilation during the 
active neutrino cooling of the accretion disk is that the magnetosphere for the 
BZ process might be disturbed so much that the BZ process is suspended until 
the burst of e+ e~ pairs clears out. However, since the pairs are produced 
with strong directionality along the rotation axis, most of the produced pairs 



expand along the axis in less than a second |2^ . This is reasonable because the 
electric field along the magnetic field lines is almost negligible and therefore 
it will take too long for the e"*" to reverse its velocity to make the same 
current as the e~. Also the strong magnetic field keeps them from moving 
perpendicular to the magnetic field lines so that the separations between 
them cannot be effective. Therefore the pair contribution to the net current 
flow can be negligible and it might not disturb the magnetosphere for the BZ 
process so violently Q. However, the effects of the neutrino cooling process 
of the accretion disk can be considered as an additional disturbing burst by 
the ui? driven e+e~ burst, which lasts less than few seconds of the BZ burst. 



^Although very hard to estimate, a very small fraction of pairs which has very small 
momentum component along the axis can contribute to BZ currents along the magnetic 
field lines, which might be sufficient for the force-free configuration. 
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The structure of a force-free niagnetosphere can be described by a stream 
function ip{r^6) and two functions of ip] VLpi^ip) and B^{ip) (or equivalently 
I = —\aujB^) 0][^. Hence ifj aX a. point of {r^O) is equal to the total 



magnetic flux upward through the azimuthal loop at (r, 0). On the other 
hand from the equation of motion, the stream function is proportional to the 
toroidal component of the vector potential, if) = 2t[A^. On the horizon, it 
determines the total magnetic flux, \E'(6'), through the horizon up to 9, ^{0) = 
ipi^TH, 0). The poloidal and toroidal components of the electromagnetic field 
are 

^P = _^k^i^V^, E^ = Q (73) 

2-na ^ ^ 

B- = Vix0^^,^_2/^ 
27100 auj 

where I{r, 9) is the total current downward through the loop at (r, 9). Qf{4') 
is the angular velocity of the magnetic field relative to absolute space. Bland- 
ford and Znajek0 derived the differential equation (stream equation) for the 
stream function, which takes the formp^ 



It includes two functions of t/^, which are not known a priori. In solving 



eq. ([75|) the boundaries of the force- free region and the boundary conditions 
on ■?/;, r2p, and / should be specified. There is no known analytical solution 
of eq. ([75|) . Blandford and Znajek obtained perturbative solutions for small 
ajM using the analytic solutions in charge free space around the non-rotating 
Schwarzschild black hole. With fip = 0,ti; = and J = 0, the stream 
equation reduces to 

V{^W} = 0. (76) 

MacDonald|2^ developed a numerical method to obtain solutions with finite 



ajM., in which the solutions of eq. (^) with appropriate boundary conditions 
are spun up numerically. It has been found that the poloidal field structure 
does not change greatly as the holes and fields are spun up. This result 
implies that the flux of the poloidal magnetic field threading black hole does 
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not decrease greatly as a ^ M in contrast to the rotating black hole in free 

space. 

One of the interesting solutions is the poloidal field structure generated by 

the paraboloidal magnetic field solution of eg. (|76D : 



i! = j7^{{r - 2) (1 - cos 0) + 2[2 In 2 - (1 + cos 6) ln(l + cos 6)] } (77) 

where ipQ is the total flux threading the hole. Since this solution extends not 
only onto the horizon but also onto the equatorial plane where the accretion 
disk is supposed to be placed, the boundary conditions on the disk should 
be satisfied by the spun up solutions generated by eq.(^). The bound- 
ary conditions depend strongly on the accretion disk model, which will be 
discussed in the next section. The presence of the accretion disk may be con- 
sidered to be the main source of the difficulty in obtaining a solution with 
these complicated boundary conditions. However, the solutions with proper 
boundary conditions provide us a way as to how we can infer the magnetic 



field from those developed in the disk. It has been demonstrated [p!0| l pT| that 
the strength of the magnetic field on the black hole is not much stronger than 
those on the inner disk. 



6 Rotation of the Black Hole 

It has been suggested that merging compact binary systems |^ or hypernovae 
[^] may result in a disk with rapidly rotating black hole at the center. In 
this section we will demonstrate how this is possible using semiquantitative 
arguments. The basic idea is that a substantial part of the orbital angu- 
lar momentum of the binary system or the spin angular momentum of the 
progenitor of hypernova can be imparted onto the black hole. 

Consider first the BH-NS merger. The typical distance for the merging 
system in this case is the tidal radius, Rt, at which the neutron star (radius 
Rns) fills the Roche lobe[p^: 



where 



Rt = :^(l±^)i/3, (78) 

0.46^ g 



^ (79) 



If the tidal radius is greater than the last stable orbit radius Q, Ris, we can 
calculate the Keplerian orbital angular velocity 



^K = ^^. (80) 

where 

M = Mbh + Mns. (81) 

Then the orbital angular momentum of the binary system can be written 

Jbinary = fiRt^K (82) 

I (^ D 

= Mbh Mns J j^ (83) 

where fx is the reduced mass, (x = MmilMms,^ Using eq. (|78D , 



J,_, = 1A7 MBHM^s^l^i^f' (84) 

During the collapse, a part of angular momentum is carried off by gravita- 
tional waves(or possibly in the later stage by neutrino cooling ) and a part of 
the total mass explodes away or remains in the torus around the black hole. 
Assuming that a fraction(a;|/) of the orbital angular momentum goes into the 
black hole, which keeps a fraction(?/) of the total mass M|, we can calculate 
the angular momentum parameter of the black hole, a: 

G 

JbH = a{yMf— = XyJunary (85) 

"" ^-^yM^^ GM ^Mns' ^ ' 

For Mbh = 2.5Mq, Mns = l-SM©, Rt = lO^cm, 

5 = 0.53- (87) 



y 



■^This leads to the condition for the black hole mass|2||: Mbh < MNs[{OA§ff£)^/'^ 
1] - 2.3Afo. 

^x is a fraction of the specific angular momentum 
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which is a quite reasonable value for an efficient Blandford-Znajek process. 

For the NS-NS merger, the radius of physical contact (2_RAr5) is smaller 
than the tidal radius, 

Rt = 2AQRns2^^^ = 3Rns- (88) 

Hence it is reasonable to consider the orbital angular momentum at the tidal 
radius. Following the same procedure as in the BH-NS merger, we get 




(89) 

where we replaced Mbh by M^s in eg . (^) . 

In the hypernova modelQI^^, a massive rapidly spinning progenitor 



{Mo ~ 4OM0) is supposed to collapse into a rapidly rotating black hole. 
If we assume a critically rotating progenitor(fi = Q^ at the surface, Ro) and 
solid body rotation throughout, the angular momentum of the core which 
eventually collapses into the black hole is 

2 

Jcore = -McoreRcore^K{Ro) (90) 

5 



where the moment of inertia of the core is assumed be that of a uniformly 

2 A/f f?2 



distributed spherical object^, / = ^McoreRcore- Using 



n. ^ ^ (.) 



we get 



7 -^M 7?2 ^^ m-\ 



As in the previous merger case, a part of the core angular momentum(x) and 
the core mass(?/) collapses into the black \io\q{Mbh = V^core)- Then we get 

JbH = XyJcore (93) 



^Since in general / = k'^McoreRcorey ^^^ ^^ << 1 fo^' radiative starspq, this is an 
upper limit to the true angular momentum. 
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xy^M^^^^Rl^^^l^ (94) 

d{yM,oref-. (95) 

c 



Then the black hole angular momentum parameter is given by 

(96) 



2x i?2 [GMo 



5yGM,ore^ Rl 



/ X , -flcore \ 3/2 l ^core^ / JV^o 



\3/2 



by' Ro \ GM,ore\ M, 



(97) 



core 



With 



M„ ~ ISMq, , M^ore ~ 2M0(~ 3fcm) (98) 



a ~ 0.1- (99) 



we get 



This depends strongly on the numbers taken in eq.(p8D and on how the core 
angular momentum is determined. The specific angular momentum of the 
core can be calculated using eg. (p^) and eq.(^) as 



2 IGM 

(^core ^ core! ^''^ core '^^coreXl p3 V^^J 

= 1.8 X lO^^cm^^s"^ (101) 

which is much smaller than that from the numerical simulation of collapsars0: 
o-core ~ lQ^^cm?s~^ . Of coursc, if the core angular momentum is not com- 
pletely redistributed during the precollapse evolution the core is likely to be 
spinning faster than VIk'-, The maximum possible value comes from replacing 
VLk determined at the progenitor radius in eq.(^2|) by the one determined at 
the core radius 



V ^core 

= (1.3 X 10^)^^, (103) 
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and we get 

a^or-e = 2.34 X 10^^cm^s-\ (104) 

Then the angular momentum parameter can be given by 

a = -^^^ = 2.3- (105) 

y ^max y 

In summary, it is very plausible to have a rapidly rotating black hole(a > 
0.1) as a resulting object in the center in the merging systems and also in 
hypernovae of large angular momentum progenitors, but a precise value of a 
will be difficult to calculate. 

7 Magnetized accretion disks 

A black hole by itself cannot keep magnetic fields on it for a long time. 
Magnetic fields diffuse away in a short time ~ i?f/-/cP]0]. The most plau- 
sible environments which can support a magnetic field threading the black 
hole are accretion disks surrounding the black hole. There are two issues 
about accretion disks that we need to consider in order to decide whether 
the Blandford-Znajek process is a viable power source for gamma ray bursts: 
The life time should be long enough to extract the bulk of the black hole 
spin energy and also the magnetic field on the disk should be strong enough 
to power the gamma ray bursts from the spinning black hole. A strong mag- 
netic field on the inner part of the accretion disk is neccessary to keep the 
magnetic pressure comparable to that of magnetic fields on the black holef^]. 
It has been also demonstrated ||11| from the axisymmetric solutions discussed 



in section 4. Since the magnetic field on the disk affects the angular momen- 



tum transfer of the accretion disk via magnetic braking ||2^ and/or magnetic 
viscosity, the accretion process is not independent of the magnetic field on 
the disk. The magnetic field on the disk should not be larger than the value 
from the equipartition argument: 

^ ~ Pa^sk, , (106) 

and Beqp can be considered as an upper limit to the magnetic field which can 
be supported by the accretion disk. Recent numerical calculations on the 
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hyper-accreting black hole by Popham et al.0 show 

Pdisk ~ lO^'^erg/cm^ (107) 

Begp ~ 5 X lO^^gauss, (108) 

which implies the accretion disk of the hyper-accreting black hole may sup- 
port a magnetic field strong enough for the gamma ray burst. However it 
depends on the detailed mechanism how strong a magnetic field can be built 
up on the disk. One of the approaches is that the magnetic field is evolved 
magnetohydrodynamically during the accreting process, which depends on 
the magnetic viscosity. 

as as 



The magnetic viscosity, u"^""^, is defined pO 



Att ^ dr ' ^ ' 

which can be parametrized by the viscosity parameter, a, in terms of dimen- 



sionful quantities of the disk asp^ 



^mag ^^^ag^^jj (110) 



where Cg is the sound velocity of the disk(cs = J^Pdisk/p) and if is a half of 
the disk thickness. Using hydrostatic equilibrium perpendicular to the disk 
plane, 



H = .^/^d.su (111) 



P 



we get 



^disk^/^ 



d 



(112) 



For a Keplerian orbit 



^mag ^ ^mag s _ ^^^^^3^ 

^disk^/7 



r^ — = --^Idisk, (114) 

dr 2 
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eq.( |109| ) can be written as 

TDdyn r)dyn 



An 



T^c^^""' Pd^skVl (115) 



We can see that the Maxwell stress of the magnetic field which has been 
built up by the accretion dynamo is also proportional to the disk pressure 
but for a different reason from that of the equipartition argument, eq. (|106|) . 



Numerical estimations of the viscosity parameter, a™"^, obtained for various 
boundary conditions range from 0.001 to 0.005|3^. This means that the 



magnetic pressure is only a small fraction of the disk pressure 0]. Using the 
estimation of disk pressure by Popham et al., we can estimate the dynamically 
generated magnetic field by the accretion: 

^dyn _ IQlS^awss (116) 

which may not be strong enough for a gamma ray bursts powered by the 
Blandford-Znajek process. 

However, an accretion dynamo might not be the only process responsible 
for the magnetic fields on the disk. We know from the recent observations 
that there are a number of pulsars, magnetars, which are believed to have 
a strong magnetic field of ~ 10^^ G. Although the origin]^ of such strong 
magnetic fields is not well known at the moment, one may consider the case 
that a debris torus or disk around the black hole which was formed by the 
disruption of a neutron star retain the high ordered field of that neutron star. 

Now consider the axisymmetric solution around a disk with a force-free 
magnetosphere which has been discussed by Blandford||29||. Here we adopt 
cylindrical coordinates, (r, 0,^;), where the 2;— direction is perpendicular to 
the disk. The sum of current fiows into the disk up to radius r defines the 
surface current density, Jr, which is proportional to the poloidal component 
of the magnetic field on the disk, 5^: 

i^ = -B, (117) 

where 2 in the denominator comes from the fact that the radial current 
density Jr includes the currents into the disk both from above and below. 
The torque exerted by the annular ring with width Ar of the disk due to the 
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Lorentz force is given by 

AT = -r2TTrJrB,Ar (118) 

= 2r^^AS,AS = 2nrAr (119) 

4:71 ^ 

For the steady state accretion disk with surface density, S, the angular mo- 
mentum conservation can be written by 

j:y^^^!^2nrAr = AT + ^Ar (120) 

or or 

M^Ar = AT+^Ar (121) 

or or 



where the torque due to the shear force of differential rotation, Gj^, is given 
by 

G = 2nruE(r^—)r (122) 

Or 

To see the effect of magnetic field, we consider only the magnetic viscosity 
for the moment. Using eg . (11091) , we get 



Qmag ^ 27Tr^-^^ (123) 

p An 
= 4^r''H^^ (124) 

= ^^r^^ML (125) 

i Idisk 47r 

where 2H p = S. Assuming a simple power dependence of H^^ oc r", 

^^ = 4(2 + n)nrH^^ (126) 

or An 

and we get from eq.( p.21| ) 

^a(r^U^ ^ ,2^ 5^[i ^ (2 + n)-^] (127) 

or r Bz 

Since H <^r, the accretion rate is determined by the magnetic braking as far 
as Br is not much larger than B^ and the second term in the right hand side 
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of the above equation can be neglected. For r ^ GMbh/c^ where the disk 
angular velocity can be approximated by a Keplerian velocity, d{r'^Q)/dr = 

M = 2rB^Bjndisk (128) 



Using the axisymmetric solution p9| 



B^ = 2rnaiskBjc (129) 

we get 

M = 4r^Bl/c (130) 

Since for a steady accretion M is independent of r, B^ ~ 1/r at large distance. 
For a numerical estimation, we take r = lO^cm and B^ = 10^'^gauss Then 
we get 

M = 6xl0-\^—^ ?MqS-^ (131) 

which is much larger than Eddington luminosity, MecI ~ 10^^^ -jf^ MqS~^. 
The total accretion during the gamma ray burst period is at most 10~^Mq, 
which may not change the discussions on Blandford-Znajek process in the 
previous sections. Hence black holes surrounded by the hyper-accretion disks 
might be a good candidate for the central engine of the gamma ray bursts. 
The magnetic field on the disk extracts not only the angular momentum as 
described above but also a substantial energy out of disk as it does on the 
black hole. The power of the disk magnetic braking can be calculated using 
eq-dHI) 

'^Pdisk = ^disk^T (132) 

= 2^^{rn,,,,)/\S (133) 

Using eq.( |129|) with Keplerian angular velocity of the disk , we get 



Bl GMbh- 



APaisk = ic^{—f^)AS (134) 



dPd^sk = ^^By^^dS, (135) 



TTC r-^ 
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then setting dS = 27irdr and using the steady-state relation B],r'^ = Bz{rin)'^rf^, 
we can obtain the total power by integrating eq. ( |135|) from Tin to infinity, 

Pdisk = 2Bz{rinfrin ^ (136) 

c 

where rj„ is the distance from the black hole to the inner edge of the accretion 
disk. Compared to the BZ power from the spinning black hole, eq . (12461) ), 
the ratio can be given by 



disk Q ' in^ / -^zy' in) \2 



f'inC / ^zy'in) s2 (^'\7\ 

Pbh ^GMfihya^^ Eh ' ^ ' 

The current conservation condition, namely that the total current flows 
onto the black hole should go into the inner edge(rj„) of the accretion disk, 
implies 

2MB^{e = 7r/2) = Cj{nn)Bf\nn). (138) 

Since the cylindrical radius in Kerr geometry u;(rj„) > 2M, we can see that 
B^ is larger than Bf^^. From the boundary conditions on the horizon |0 in 
the optimal case, 

B^ = -nnMBn, (139) 



and on the accretion diskp9| with angular velocity fi^), 

^disk ^ _2fi^r5f ^^ (140) 

we get 

B.(.rt„} [GMaGM 

and the power of the disk magnetic braking Pdisk = -jf^i^^YPBZ- It shows 
that the magnetic field on the horizon cannot be smaller than that on the 



inner edge of the accretion disk[^ and the disk power need not be substan- 



tially larger than that from the black hole[|TOl UTI I. 



8 Conclusion 

We have evaluated the power and energy that can be extracted from a rotat- 
ing black hole immersed in a magnetic field, the Blandford-Znajek effect. We 
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improve on earlier calculations to find that the power from a black hole of 
given mass immersed a an external field is ten times greater than previously 
thought. The amount of energy that can be extracted from a black hole in 
this way is limited by the fact that only 29% of the rest mass of a black hole 
can be in rotational energy, and that the maximum efficiency with which 
energy can be extracted from the hole via the Blanford-Znajek effect is 31%. 
The net amount of energy that can be extracted is therefore 9% of the rest 
energy of the black hole. We consider various scenarios for the formation of 
rotating black holes in gamma-ray burst engines, and while the resulting an- 
gular momenta are quite uncertain in some cases, it seems that the required 
values of the rotation parameter, a ^ 0.5, are achievable. 

The rate at which angular momentum is extracted depends on the mag- 
netic field applied to the hole. A field of 10^^ G will extract the energy in 
less than 1000 s, so time scales typical of gamma-ray bursts can be obtained. 
Since the black hole cannot carry a field, there must be an ambient gas in 
which the field is anchored that drives the Poynting fiux from the black hole. 
The most obvious place for it is the accretion disk or debris torus surrounding 
the black hole just after it formed. 

It has been argued that a field turbulently generated in the disk would not 
give a strong Blandford-Znajek fiux||TO[|, because the disk would dominate the 



total Poynting output, and no more than the disk's binding energy could be 
extracted. We show explicitly that a field in the disk could be much greater, 
for example if it is derived from the large, ordered field of a neutron star that 



was disrupted. The field distribution proposed by Blandfordp9| for such a 
case would allow the field on the hole to be much greater than on the disk, 
such that the Poynting fiow would not be dominated by the disk and not 
subject to any obvious limits imposed by the disk. 

We also note that a Poynting fiow may provide an alternative way of 
providing a very large magnetic field for the shocked material that radiates 
the afterglow and the gamma-ray burst itself: the standard assumption is 
that the required high fields grow turbulently in the shocked gas, up to near- 
equipartition values. But the field in the Poynting fiow only decreases as 1/r, 
so if it is 10^^ G at r = 10^ cm, it could be as high as 10^ G at the deceleration 
radius (10^^ cm), ample to cause an energetic gamma-ray burst. 
Acknowledgement 

We would like to thank Roger Blandford, Sterl Phinney and Kip Thorne 
for guidance and useful discussions. HKL is very grateful to the Nuclear 

28 



Theory Group, SUNY at Stony Brook for the hospitahty during his one year 
stay, where the most of collaborated works on this paper has been done.. This 
work is supported by the U.S. Department of Energy under Grant No. DE- 
FG02-88ER40388. HKL is also supported in part by KOSEF-985-0200-001-2 
and BSRI-98-2441. 



29 



Appendix 1. Rotating Black Hole 

A rotating black hole is defined as its angular monientuni( J) and niass(M). 
The angular momentum is measured by the non-newtonian gravitational ef- 
fect (gravitomagnetic effect) far from the black hole and the mass is mea- 
sured by the newtonian gravitational field far from the black hole. The solu- 



tion of Einstein's equationl^S] defines the metric around the rotating black 
hole with specific angular momentum a = J/M. Using the Boyer-Lindquist 
coordinates [0, in the natural unit G = c = 1 the metric can be written as 

a and /3 are lapse function and shift function defined respectively as 

a = ^ (143) 

^* - -^ (144) 

/?^ = /3" = (145) 

(146) 



where 



A = 


r^ + 0^ - 2Mr 


P' = 


r^ + a^ cos^ 9 


Y? = 


(P + a^f-a'Asm'e 


The metric tensor gij is given 


by 


9rr 


P^ 2 ~2 

-^, 9ee — P , 94,4) — ^ 


where 






Y, 
uj = — sin 6^ 
P 



(147) 
(148) 
(149) 
(150) 



(151) 



(152) 



The volume and area elements dV and dS are defined in the standard wayp 
by 



dV = J det{gij)drd6d(j) (153) 
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2 ~ 

^ ^-.drded(j) (154) 



I Ox"^ Ox 

dSi = \Jdet{gij)eijk-Kr-K-dbdc (155) 

The circumference of a circle around the axis of symmetry is 2TT,/g^ = 2ttuj. 
The horizon, rn, is defined as a larger root of Ar=r„ = 0, 

rjj = M+{M'^-aY^^ (156) 

2MrH = rl + a^ (157) 

The surface area of the horizon, Ah, is given by 

Ah = 4n{rl + a^) (158) 

and the entropy and the irreducible mass of the black hole are given respec- 
tively by 

Sh = ^Ah (159) 

= -^-{^H + a) = -——Mth- (160) 

n n 



M.„ = ^^=i,^r7;i=^|£ (161) 

The the mass of the black hole can be rewritten in terms of J, Mirr, and Sh, 



The angular velocity of the black hole is given by 



J 1 



a2 



2M3(1 + ^1- JVM4) 2Ml + v^ 
where a is the angular momentum parameter defined as a = J/M"^ . 



(164) 
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Appendix 2. Electromagnetic fields in vaccum arround rotating 
black hole 

For a rotating black hole immersed in an asymptotically uniform magnetic 
field 

B = Bz, r^oo, (165) 

the poloidal components of electric and magnetic fields can be written ana- 
lytically § 

B dX , , 

B^ dX , , 

«» = -2E^* («^' 

-BaE.da^ Msin^^,. ,2..^^,r. ,^^„, 

„ —BaJ^.da^ sin^6',^o . 9 , , x c? . 1 ^ /,^„n 

where X = 2iS^(S^ — Aa^Mr) and a is the specific angular momentum of 
the black hole, J = aM. Here we adopt the natural units G = c = 1. 

The radial components of magnetic field on the horizon, B^ , can be 
written explicitly by 

' " (r|, + (a/c)2cos2^)2 ^ ' 

\{r\ - {alcf){r\ - {ajcf cos' d) + 2(a/c) V^lrj^ - M)(l + cos' Q)\ 

T^H ^ ^^ n7T) 

2r^c(r|, + (a/c)2 cos2 ^)2 ^ ' 

{rl - {ajcf) [(3r| + (a/c)'(l + cos' ^)) cos' d - r\\ 

(172) 

It is clear that B^ vanishes as the rotation of black hole approaches extreme 
rotation 

a^GMjc, th^GM/c^ (173) 



This is what has been observed [^ pO| as the absence of the magnetic fiux 



across the maximally rotating black hole. There are no (j) components of 
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electric and magnetic fields in eq.( |169| ), and the 6 components vanish as v A 
as the horizon is approached. The main reason for these vanishing tangential 
components is that there are no currents outside the horizon in charge-free 
space. In the absence of in/out currents there is no toroidal magnetic field on 
the horizon. This also requires a vanishing tangential component of electric 
field on the horizon, which induces a charge separation on the horizon]^]. 
Therefore there is no current on which the external magnetic field exerts 
torques to slow the black hole down so as to extract its rotational energy. 
One can verify from eq.( [166| ) - eq. (|169|) that there is no outward component 



of the Poynting vector (essentially it is due to the absence of the component 
of the magnetic field) , which means no outflow of energy to infinity. 
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Appendix 3. Axial-symmetric Force-free Magnetosphere 

In this appendix, the original forniulation|0 of the Blandford-Znajek pro- 
cess is summerized. Consider the stationary and axial symmetric situation 
where all the partial derivatives of the physical quantities with respect to 
time and azimuthal angle are vanishing. The electromagnetic field tensor 
is given by the vector potential A^ as 

F^u = A,,^-A^,. (174) 

(175) 

and the magnetic field is given by 

(176) 
(177) 
(178) 

The magnetic fiux, \E', through a circuit encircling = — >■ 27r, 

^ = JA^d(t) = 2TiA^ (179) 

It defines a magnetic surface on which A^{r, 9) is constant and therefore is 

characterized by the magnetic fiux \I' contained inside it. Magnetic field lines 

are spiraling on the surface and no magnetic field lines can cross the magnetic 

surface. 

(1) Force- free magnetosphere 

The force-free condition for the magnetosphere with the current density 



Br 


= 


Fetf, = 


= A^fi 




Be 


= 


F^r = 


- ^0,r 




Bs 


= 


Fre = 


= Ae^r - 


A; 



can be written by 



F^^J^ = 0, (180) 



Ao,,, J'- + v4o,e J^ = (181) 

A^^rr + A^^eJ' = (182) 

-Ao,rJ^-A^^rJ''' = B^J'^ (183) 

A^^eJ"" + A^^eJ'*' = B^r. (184) 
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/,From eq.(^) and ( [I82D we get 



(185) 



We can see that Aq is also constant along the magnetic field lines and the 
electric field is always perpendicular to the magnetic surface. We can define 
a function u^r, 9), 

= =-u (186) 

dAo = -oodA^, (187) 

which is also constant along the magnetic surface, uj can be identified as an 
electromagnetic angular velocity(an angular velocity of magnetic field line on 
a magnetic surface). We can also suppose a function fi{r,9) satisfying 

An.r = -^A^,e 

A-nf = —Asr (188) 



where ^/g = J —det{g^y) = p sm6. The current conservation 

-)^(v^^"),. = (189) 

leads to 

{nA^,e),r = {pA^,r),e (190) 

which implies that /i is also constant on the magnetic surface. 
From eq.( p.83| )(or equivalently from eq.( |184| )), we get 



j^ = ^jo + -J^B^ (191) 



Together with eq. (|188| ) , we now have an expression of the current in terms of 
the magnetic field and charge density with yet undetermined function u and 
/i. The outward current / can be calculated from the r and 6 components of 

35 



the current. Since the outward curent is proportinal to the magnetic field, 
the current between the magnetic surface, dJ, can be written in terms of the 
magnetic fiux as 

dl = —^id^i (192) 

= -^^dA^ (193) 

(2) Energy and Angular Momentum Outflow 

For a stationry and axial symmetric system, the conserved fiuxes can be 
defined about the axis of symmetry. The force-free condition ensures the 
conserved electromagnetic energy fiux 

^^ = T^ (194) 

and angular momentum fiux 

£^ = -T^ (195) 

where the electromagnetic energy momentum tensor is given by 






T^^ = —[F^^F^P - -g^'Tp^pP''). (196) 



The energy fiux and angular momentum fiux have a simple relation, 

S^" = cj£^ (197) 

and the poloidal components can be written as 

£" = -^^^A^,eB^ (198) 

47r p^ 

I A 

UJ — 

An p 
Using the physical component, £f 



£' = —u-A^^rB^ (199) 



£-f — \J gm/grr£-r — \^9oo9^£^ (200) 
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we get the power at infinity given by 

1 



9oo 

4n I '^\l~ 



1 

in 



UJa 



grrQee 



Bj,Aa 



~TI / ^( 2 — B^)A^^0d6d(j) 



Defining -Bt[H as 



B. 



A sin 9 



B, 



the power can be written as 

P = —juoBrA^^ede 



1 



ujExdA^ 



From the inhoniogeneous Maxwell equations 



we get 



Using eq. (|188|) 



Bxfi — 


An^r 


BT,r = 


-Ati^J 


fJ'A^^e 


= Bxfi 


f^A^^r 


= Bt^t 



(201) 
(202) 

(203) 
(204) 

(205) 

(206) 
(207) 

(208) 



(209) 
(210) 



(211) 
which shows that Bt is also constant along the magnetic surface. Since 

lidA^ = dBr, (212) 
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the outward current can be calculated as 



and we get 



dl = ^fxdA^ (213) 

= ^dBr (214) 



m = \jdBT=\BT{e) (215) 



where Bt{0 = 0) = has been used. Using the physical component , -B; 



B. 







eg. ( pi 51 ) can be written 

''-fiaBi (216) 



y/grrgee '■ 

which is nothing but the Ampere's law on the stretched horizon 

2TiuBf = 47r/ (217) 

using the tangential field on the stretched horizon]^ 

Bf = aB^ (218) 

Using eq. (^), (|?TB| ) and (^) we get 

dP = -cuIdA^ (219) 

= cu/rf^ (220) 

2n 

P = -— I uldm (221) 

27r J 

Since w, / and ^ are constant along the magnetic surface we can evaluate 
the integral on the horizon of the black hole: 

P = / ojld-^ (222) 

27r J horizon 

Hence we can see that it is a proof of the simple-minded derivation of eq.(21) 
in the text with a; identified with Qp. 
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Appendix 4. Power from Rotating Black Hole 

Using the boundary condition ^^ tliat Br is finite and 



sin6'[u;(r|^ + a 



2^ 



a\ 



Ph 
which can be written for B^ as 



Bt = ^^vS -^Br (223) 

rfj + a'' cos 

- -M^tiZ^EEHB^ (224) 



A p 



B, = L_ y---^^>-^ Br (225) 



A 



Br (226) 



^"-^^)%^5, (227) 



A 
and the toroidal component on the stretched horizon is 

5f = (a; - nH)^BH (228) 



where Bh = Bf. Then the current up to 6* = 9h can be written using eq.(pT 

I{9H) = ~sm9B^ = ^{u-nH)u^BH (229) 

which is exactly the same current derived in the membrane paradigm|P]. 

The total power out of the black hole can be calculated using eq. ( p06| ), 
(|T5D, and (|2^) , 

P = - I uBrBrdO (230) 

= - j ujIpuBude (231) 

= -- j uj{uj~VLH)B]jpCj^de (232) 

Adopting the optimal condition0, uj ~ ^^72, the power can be written as 

P = -- fuj{uj-nH)B%pCu^d9 (233) 
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r% + a^Ysin^ 9 



-^^2M^^ I rl + a-cos-9 ^^^^ ^^^^^ 



32'M' r|^ ' ' a ' '7 (^)2 + cos2^ 
Using the integral identity, 

where /i is defined in ref [H^ as 

/i = — (237) 

J/M^ " 



1 + y^l - jyM^ l + VT 
M _ M a _h 
th a th a 



a2- 



(238) 



(239) 



we integrate the angular integral assuming the average magnetic field, < 



1 a ^ {2MrH?{l + h'') 
32 W^ rl h^ 



< Bjj>2[{h + 1/h) arctan/i - 1^40) 



H 

\{^?M'<Bl>f{h) (241) 



where 



f{h) = ,0^ M (/i + l//i)arctan/i-l] (242) 

-^ 2/3, for h^O, (243) 

= vr — 2, /or h = 1( extreme rotation, a = M) (244) 

Now inserting G and c properly, we get 

^ = 3<^''^'' < "" > ^<">|r (245) 

= hHl-'<Bj,>f(h)^ (246) 

4 c"* 
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Numerically, 

^MUlO^'gaussf = (2 x 10^^ x lO^^f x i^^^^^^ (248) 

= 6.7 X 10^° er^/s (249) 

Then the total power of Poynting flux is 

P.l.TxlO»(^«^«^^)V(/.)e.«/. (250, 

One of the frequently-quoted numbers for the total power of BZ process 
from ref . [Q is 

which is one order of magnitude lower than eq. ( ^45| ) Q. In ref [^ , the angular 



integration of eq. (|222|) is approximated by replacing 





u- -. rl/2 


(253) 




A^ -> -ktIBh 


(254) 


p 


- usC/^''<'''''>^'^^'"€ 


(255) 


g\ac/GM) 


= l + (l-(ac/GM)2)V2 


(256) 




-^ 1 for a ^ 


(257) 




= 2, for a = M{ extreme rotation) 


(258) 



Then 



P(eg.255) g^{ac/GM) 



P{eq.^) 32f{h) ^^^^^ 

-^ 3/64 Jora -^ (260) 

~ 5/56 , for a = M{ extreme rotation) (261) 

which explains the difference of an order-of-magnitude. 



Recently Popham et al. M use 



P-^0''i^)'i lf"^ )'^rg/s (252) 



for M = 3Mq hence equivalent to eq.(E51J). 
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